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Abstract. We analyze the properties of the energy landscape of finite-size fully 
connected p-spin-like models. In the thermodynamic limit the high temperature 
phase is described by the schematic Mode Coupling Theory of super-cooled liquids. 
In this limit the barriers between different basins are infinite below the critical 
dynamical temperature the ergodicity is broken on in infinite times. We show 
that finite-size fully connected p-spin-like models, where activated processes are 
possible, do exhibit properties typical of real super-cooled liquid when both are 
near the critical glass transition. Our results support the conclusion that fully- 
connected p-spin-like models are the natural statistical mechanical models for 
studying the glass transition in super-cooled liquids. 
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Introduction 

In recent years a significant effort has been devoted to understand slow relaxation 
dynamics observed in many, apparently unrelated, systems such as structural glasses, 
spin glasses, disordered and granular materials or proteins among others. In such 
systems the characteristic relaxation time may change of many orders of magnitude if 
the external parameters, e.g. the temperature T, are slightly varied. As a consequence 
correlations display non-exponential behavior, and equilibration processes slow down 
giving rise to non-equilibrium phenomena known as aging. 

The common denominator which makes all these systems displaying similar 
behavior near the (dynamical) critical temperature is the complexity of the energy 
landscape. The trajectory of the representative point in the configuration space can 
be viewed as a path in a multidimensional potential energy surface . The dynamics 
is therefore strongly influenced by the topography of the potential energy landscape: 
local minima, barriers, basins of attraction an other topological properties all influence 
the dynamics. 

The potential energy surface of a super-cooled liquid contains a large number of 
local minima, called inherent structures (IS) by Stillinger g]. All states that under 
local energy minimization will flow into the same IS define the basin of the IS (valley) . 
With this pictures in mind the time evolution of the system can be seen as the result 
of two different processes: thermal relaxation into basins [intra-basin motion) and 
thermally activated potential energy barrier crossing between different basins (inter- 
basin motion). When the temperature is lowered down to the order of the critical 
Mode Coupling Theory (MCT) temperature Tmct the inter-basin motion slows down 
and the relaxation dynamics is dominated by the slow thermally activated crossing of 
potential energy barriers || . If the temperature is further reduced the relaxation time 
eventually becomes of the same order of the observation time and the system falls out 
of equilibrium since there is not enough time to cross barriers and equilibrate. This 
define the "experimental" glass transition temperature T g . 

The regime between Tmct and T g cannot be described by the MCT since it 
neglects activated processes responsible for barrier crossing. In MCT the relaxation 
time diverges at Tmct-, leading to T g = Tmct, and the dynamics remains confined 
into a single basin forever. Attempts to overcome these difficulties in MCT have been 
done, but probably the most clear picture comes from some spin-glass models. The 
essential features of MCT for glass-forming systems are also common to some fully 
connected spin glass models [|J, the most well known being the spherical p-spin spin 
glass model [[| We shall call these models mean-field p- spin-like glass models. The 
central point is that near Tmct the behavior of the system is mainly due to the IS 
organization (density, basins, barriers and so on), so that all systems with similar IS 
structure should have similar critical behavior. 

In the thermodynamic limit the high-temperature phase - paramagnetic in the 
spin-glass language and liquid in glass language - of mean-field p-spin-like is described 
by the schematic MCT for super-cooled liquids |t], || . As a consequence at the critical 
temperature, called Tjj in p-spin language, an ergodic to non-ergodic transition takes 
place. Below this temperature the system is dynamically confined to a metastable state 
(a basin) J8| since relaxation to true equilibrium can only take place via activated 
processes, absent in mean-field models. Therefore in mean-field models, similar to 
what happens in MCT, at To the relaxation time diverges. For these systems, 
nevertheless, it is known that the true equilibrium transition to a low temperature 
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phase occurs below Tjj at the static critical temperature T c , also denoted by T\ rs ^ 
||. This is the analogous of the Kauzmann temperature for liquids. The glass 
transition temperature T g of real systems sits somewhere in between T c and Tjj . This 
transition, obviously, cannot be reached even on infinite time in mean-field models. 

Despite these difficulties mean-field models, having the clear advantage of being 
analytically tractable, have been largely used to study the properties of fragile 
glassy systems, especially between the dynamical temperature Try and the static 
temperature T c where a real thermodynamic phase transition driven by the collapse 
of the configurational entropy takes place. The picture that emerges is however not 
complete since activated process cannot be captured by mean-field models. Therefore 
the relevance of mean-field results for real systems it is still not completely stated. 
Let us also remark that, despite the large amount of analytical work devoted to the 
study of the static as well as dynamical properties in the N — > oo limit, much less is 
known concerning the finite N behavior. 

In this work we investigate numerically finite-size fully-connected p-spin-like 
models, where activated processes are present. Comparing our results with the 
observed behavior of super-cooled liquids near Tmct we conclude that, once activated 
process are allowed, mean-field p-spin-like models are highly valuable for a deep 
understanding of the glass transition in real systems. 

All results reported here are for the Ising-spin Random Orthogonal Model (ROM) 
JTo[ |i"l| , however similar results are obtained using other p-spin-like models, as for 
example the Bernasconi Model and the Ising p-spin model. The advantage of ROM 
lies in its interaction term which is two-body, at difference with the p-body interaction 
of p-spin models, reducing computer memory problems. Moreover the use of Ising 
spins instead of continuous spins, as for the spherical p-spin model, allows for a larger 
configuration entropy and faster algorithms. Preliminary results for the spherical p- 
spin model are in qualitative agreement with those reported here. 



1. Thermodynamics of Inherent Structures: How to evaluate the 
configurational entropy 

The Random Orthogonal Model model is defined by the Hamiltonian (r], 11 , 

H = -2Y / J ij a i a j (1) 

ij 

where <7j = ±1 are N Ising spin variables, and is a N x N random symmetric 
orthogonal matrix with Ju = 0. Numerical simulations are performed using the Monte 
Carlo method with the Glauber algorithm. For N — > oo this model has the same 
thermodynamic properties of the p-spin model: a dynamical transition at To = 0.536, 
with threshold energy per spin eth — Eth/N — —1.87, and a static transition at 
T c = 0.25, with critical energy per spin ei rs b = —1.936 |r], [il]| . 

The TAP analysis || |ll| reveals that the phase space visited is composed by an 
exponentially large (in N) number of different basins, each labelled by the value of the 
energy density e at T = 0, separated by infinitely large (for N — > oo) barriers. The 
free energy of the e-TAP solution describes the thermodynamics within the basins 
labelled by e and at T = coincides with the local minima potential energy, i.e., the 
IS of the system. The dynamical transition is associated with IS having the largest 
basin of attraction for N — > oo, while the static transition with IS having the lowest 
accessible free energy || |l2| . 
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In the mean- field limit, the allowed values of e are between ei rs f, and e t h- Solutions 
with e larger than e t h are unstable (saddles), while solutions with e smaller than ei rs {, 
have negligible statistical weight. Moreover in the N — > oo limit IS with e = e t h 
attract most (exponentially in N) of the states and dominate the behavior of the 
system. Other IS are irrelevant for N — * oo. 

For finite N the scenario is different since not only the basins of IS with e < eth 
acquire statistical weight, but it may happen that solutions with e > eth and few 
negative directions (saddles with few downhill directions) become stable, simply 
because there are not enough degrees of freedom to hit them. 
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Figure 1. Temperature dependence of (e(T)} for N = 48 (square), N = 300 
(circle) and N = 1000 (triangle). The average is over 10 3 different equilibrium 
configurations at temperature T. The horizontal line is the N — -> oo limit. The 
arrows indicate the critical temperatures and T c (see text). The dotted line 
is the curve obtained from the configurational entropy for large N. Inset: size 
dependence of (e(T = 3)) as a function N~ a with a = 0.2 extrapolated down to 
the iV — > oo theoretical result —1.87 (triangle). 



To get more insight the IS-structure of finite systems we follow Stillinger and 
Weber [jl3| and decompose the partition sum into a sum over basins of different IS 
and a sum within each basin. Collecting all IS with the same energy e, denoting with 
exp[iVs c (_E)] de the number of IS with energy between e and e + de, and shifting the 
en ergy of each basin with that of the associated IS, the partition sum can be rewritten 

Z N (T) ~ / de exp N [-0e + s c (e) - 0/(0, e)] (2) 

where /(/?, e) can be seen as the free energy density of the system when confined in one 
of the basin associated with IS of energy e. The function s c (e) is the configurational 
entropy density also called complexity. From eq. (|^) we easily obtain the probability 
that an equilibrium configuration at temperature T = 1//3 lies in a basin associated 
with IS of energy between e and e + de: 

P N (e, T) - exp N [-0e + s c (e) - /?/(/?, e)] /Z N (T). (3) 
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Taking the N — > oo limit we recover the mean- field results 

^From the partition function (|J) we can easily compute the average internal energy 
density given by u{T) = (e + d{(3f)/d(3) = (e(T)) + (Ae(T)). The first term is the 
average energy of the IS relevant for the thermodynamics at temperature T, while 
the second term is the contribution from fluctuations inside the basin of the IS. The 
average is taken with the weight (|J) . 



0.3 




Figure 2. (a) Configurational entropy as a function of energy. The data are from 
system sizes N = 48 (empty circle) and N = 300 (filled circle), and temperatures 
T = 0.4, 0.5, 0.6, 0.7, 0.8, 0.9 and 1.0. For each curve the unknown constant has 
bee n fixed to maximize the overlap between the data and the theoretical result 
p.]] . The line is the quadratic best-fit. (b) Configurational entropy density as a 
function of temperature. The line is the result from the best-fit of s<=(e) while the 
symbols are the results from the temperature integration of eq. (ph for N = 48 
(empty circle), N = 300 (empty triangle) and N = 1000 (filled circle). 

Since we are interested into the IS structure we shall concentrate on (e(T)}. In 
the limit iV — > oo the only relevant IS are those with e = eth, and lim7v^oo(e(T)) = &th 
for any T > To- To measure (e(T)) for finite N we perform the following experiment. 
First we equilibrate the system at a given temperature T, then starting from an 
equilibrium configuration we instantaneously quench it down to T = 0. This is 
obtained by decreasing the energy along the steepest descent path. In this way we can 
identify the energy of the IS visited by the equilibrium trajectory. The experiment 
is repeated several times starting from uncorrelated equilibrium configurations at T 
and the average IS energy is computed. In figure |l| we report (e(T)) as a function of 
temperature T for system sizes N — 48, 300 and 1000. As expected, as N increases 
(e(T)) tends towards e t h- From the numerical data we found that the plateau energy, 
approaches eth with the power law (e p i a t e au) — eth ~ A^ 2 , see inset Fig. |l| Note 
that since N is finite we can equilibrate the system also below Td , down to the glassy 
transition T g (N), about 0.35 for TV = 48 and 0.5 for TV = 300, below which the system 
falls out of equilibrium jlij . 
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The figure shows that for finite N and T not too close to Tjj the thermodynamics is 
dominated by IS with e > e t h- This is more evident from the (equilibrium) probability 
distribution of e since it is centered about (e(T)) indicating that IS with e ~ (e(T)) 
have the largest basins. This scenario has been also observed in real glass- forming 
systems jT| |^ 0, [TsJ, fl9). 

^From the knowledge of IS-energy distribution we can reconstruct the complexity 
s c (e). In the temperature range where eq. (|^) is valid, we have 

fl c (e) =hxP N (e,T) +/3e +0ftf, e)+\nZ N (T) (4) 

If energy dependence of /(/?, e) can be neglected, then it is possible to superimpose 
the curves for different temperatures. The resulting curve is, except for an unknown 
constant, the complexity s c (e). The curves obtained for system sizes A" = 48 and 
300 and various temperatures between T = 0.4 and T — 1.0 are shown in figure || 
(a). The data collapse is rather good for e < —1.8. Above the curves cannot be 
superimposed anymore indicating that the energy dependence of /(/?, e) cannot be 
neglected. In liquid this is called the anharmonic threshold ||(], Q . To compare the 
result with the known analytical predictions for the ROM each curve in the figure has 
been translated to maximize the overlap with the theoretical prediction for s c (e) jnj. 
The dotted line is the quadratic best fit from which we can estimate the critical energy 
e c ~ —1.944 as the value where s c (e) vanishes in good agreement with the theoretical 
result ei rs b = —1.936 JDJ. 

Direct consequence of /(/?, e) ~ f((3) for e < —1.8 is that in this range the 
partition function can be written as the product of an intra-basin ]2^ | contribution 
[exp(— Nf3f)\ and of a configurational contribution which depends only on the IS 
energy densities distribution. The system can then be considered as composed by two 
independent subsystems: the intra-basin subsystem describing the equilibrium when 
confined within basins, and the IS subsystem describing equilibrium via activated 
processes between different basins. As the temperature is lowered and/or A" increased 
the two processes get more separated in time and the separation into two subsystems 
becomes more and more accurate. A scenario typical of super-cooled liquids near the 
MCT transition || §. 

The form of /(/?) for the specific system can be computed studying the motion 
near the IS, for example using an harmonic approximation (2^, [^J. However, this 
usually gives only a small corrections to thermodynamic quantities for T close to Tu 
fizdj , |lj and we do not consider it here. 

Another important consequence of the separation into two subsystems is that in 
ec l- © fifl\ e ) can be neglected since it cancels with the equal term coming from the 
denominator. Therefore from the knowledge of s c (e) we can easily compute the average 
IS energy density (e(T)). For large N this is given by the saddle point estimate, see 
eq. (|: 

max [—/3 e + s c (e)} . (5) 

e 

The result obtained using for s c (e) the quadratic best-fit of figure || (a) is the dotted 
line shown in Fig. |l[ The agreement with the direct numerical data is good already 
for N = 300. ^From the form of (e(T)) for large A^ we can identify the static critical 
temperature T c ~ 0.3 as the temperature below which (e(T)) remains constant, not 
far from the theoretical results T c = 0.25 plf . 

To have a check of our results we have computed the configurational entropy 
density using a different approach When the IS subsystem is in thermal 
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equilibrium then the temperature dependence of the configurational entropy density 
can be evaluated from the thermodynamic relation 

M = i (6) 
d(e(T)) T y ' 

integrating the T dependence of d{e(T))/T. Using the data of figure [l] we obtain 
the curves shown in figure || (b). The line is the result valid for large N obtained 
from the quadratic best fit of s c (e) [figure || (a)]. The agreement for N = 300 and 
1000 is rather good. Note that increasing N reduces the IS energy range explored by 
the system for a given fixed Montecarlo simulation length. This in turn reduces the 
temperature range in figure ^ (b) . 

To have a statistical description of the energy landscape we also investigated cross 
correlations among IS and equilibrium configurations at T. It emerges that while the 
equilibrium configuration is highly correlated with the corresponding IS, different IS 
are uncorrelated. Moreover the probability distribution of IS-IS overlaps is a Gaussian 
centered in zero and variance which goes to zero as N increases. The analysis of triplet 
of IS does not reveal any particular organization of states. This result, characteristic 
of the Random Energy Model is known to hold also for multispin interaction 
spin-glass models. 
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Figure 3. Average inherent structure energy in equilibrium as a function of 
temperature (a) and as a function of time during the non-equilibrium process (b). 
The system size is N = 300, Tj = 3.0 and (top to bottom) T f = 0.1, 0.2, 0.3 and 
0.4 (panel b). 



2. Non-equilibrium behavior: the role of activated processes 

More informations on the IS structure can be obtained from non-equilibrium relaxation 
processes. To study the non-equilibrium dynamics we quench at time zero the system 
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from an initial equilibrium configuration at temperature T, > T g to a final temperature 
Tf < T g and study the evolution of the average IS energy per spin (e(t)) as function of 
time. This is obtained by regularly quenching the system down to T = to calculate 
the instantaneous IS energy. The result for a system of N = 300 spins and initial 
temperature Tj = 3 is shown in figure ||] (b) for final temperature Tf = 0.1, 0.2, 0.3 
and 0.4. The average is over different equilibrium initial configurations at T. 

The analysis of the figure reveals that the relaxation process can be divided 
into two different regimes. A first regime independent of Tf, and a second regime 
independent of both Ti and Tf. The final temperature Tf controls the cross-over 
between the two regimes. A similar behavior has been observed in molecular dynamics 
simulations of super-cooled liquids Q . We note that in our case since we use discrete 
variables, and hence a faster dynamics, we do not have the very-early regime observed 
in (l^] where (e(t)) is almost independent of t. 

The two regimes are associated with different relaxation processes. In the first 
part the system has enough energy and relaxation is mainly due to path search out 
of basins through saddles of energy lower than ksTf. This part depends only on the 
initial equilibrium temperature Tj since it sets the initial phase space region. Different 
T leads to different power law. In particular relaxation must slow down as Tj decreases 
since we expect that lower states are surrounded by higher barriers. This expectation 
is supported by our numerical data. In figure ||| we report the behavior of (e(t)) as 
function of time for different initial temperatures. The slowing own of the first regime 
is clearly seen. 




Figure 4. Average inherent structure energy as function of time for initial 
temperatures Tj = 3.0, 2.0 and 0.9, final temperatures TV = 0.1, 0.2 and 0.3. 
The average is over 300 initial configurations. The system size is N = 300. The 
line denotes the slope —0.0025. 



During this process the system explores deeper and deeper valleys (basins) while 
decreasing its energy. The process stops when all barrier heights become of 0{ksTf). 
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/From now on the relaxation proceeds only via activated process. A first consequence 
is that lower the final temperature Tf shorter the first relaxation, in agreement with 
our findings [See figures |and|. 

The analysis of the distance between the instantaneous system state and the 
corresponding IS, counting the number of single spin flip needed to reach the IS, reveals 
that for all times the systems stays in configurations few spin flips away from an IS. 
The number ranging from 8 — 9 for short times to 1 — 2 at larger times. A similar 
study starting from equilibrium configurations at temperature T e ((e(i))) evaluated 
comparing panels (a) and (b) of figure || jl9) leads to similar numbers. We then 
conclude that during relaxation the aging system explores the same type of minima 
(and basins) visited in equilibrium at temperature T e . Direct consequence is that once 
the system has reached the activated regime there cannot be memory of the initial Tj , 
and all curves with different Tj but same Tf should collapse for large time [figure . 

To have more confidence with this picture of relaxation we have studied the 
distance from the instantaneous configurations and the nearest saddle to a different 
basin. This is done by counting the number of spin flips needed to reach the saddle 
through the less steep path, i.e., the path that gives the minimum increase of energy 
at each spin flip. A typical result is shown in figure ^. Note the strong slowing down 
at the "kink" where the relaxation law changes. 




Figure 5. (a) Number of single spin flip to reach the nearest saddle from the IS 
as a function of tt. (b) (e(t)) as a function of tduring the non-equilibrium process. 
The system size is N = 300, T t = 3.0 and Tf = 0.2. 

Finally in figure ^| it is shown the average IS energy of IS reached by crossing 
the saddle nearest to the instantaneous configuration. Note that for short time the 
crossing leads to IS with similar energies. By comparison we also report the analogous 
curve obtained starting from the instantaneous IS. In this case the energy is much 
lower. The picture changes near the kink where both curves merges. Moreover for 
longer times both curves are above the true relaxation curve, indicating that the 
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system does not relaxes passing through the less steep path. Indeed this is a very 
special path which may be difficult to find. Higher saddles require higher activation 
energy but can be found more easily and dominate the relaxation dynamics. More 
detailed studies of barriers heights and exit time are in progress. 




Figure 6. IS energy as function of t for N = 300, Ti = 3.0 and T f = 0.2. The 
average is over 300 initial configurations. Symbols are: circle: IS energy from 
instantaneous configuration; triangle down: IS energy crossing nearest saddle to 
instantaneous configuration; triangle up: IS energy crossing nearest saddle to 
instantaneous IS. 



3. Conclusions 



To summarize, in this work we have shown that finite-size mean-field p-spin-likc 
models are good candidates for studying the glass transition. The key point is 
that near the glass transition the thermodynamics of the systems is dominated by 
the IS distributions, therefore all systems with similar IS distributions should have 
similar behavior. Finite-size mean-field p-spin-like models have the double advantage 
of being analytically tractable for N — > oo and easily simulated numerically for finite 
N, offering good models to analyze the glass transition. Finally we note that the 
analysis presented in this paper opens the way for the study of generic glass and 
spin-glass models (such as the Edwards- Anderson model) where a careful study of the 
thermodynamics associated to the IS has never been considered. 
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